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1. Introduction 

We study a hadronic component of extensive air showers. Key parameters that we want to 
determine are a shape of a muonic subshower profile and an atmospheric depth of its maximum 
(^max)- The Xmax is the well known observable that can be used for a mass determination of 
primary cosmic ray particles, see e.g. Ref. ||I|] where detailed description of muon production 
depths reconstruction at the Pierre Auger Observatory together with mass composition implications 
of the measurement are summarized. 

There are two methods that can be used for calculation of shower profiles. The most precise of 
them is the use of Monte Carlo simulation tools that provide detailed information about showers. 
A disadvantage of using Monte Carlo tools lays in the fact that they are computationally very 
demanding. On the other hand, for a qualitative description and better physical insight one could 
use simple analytical or semi-analytical models that are in general easy to handle and have a clear 
interpretation. We focus on the latter option. 

The best known analytical model for an electromagnetic component of cosmic ray showers is 
the Heitler model [^. The Heitler model can be extended to the Heitler-Matthews model [||] to 
incorporate hadronic (pionic) component. The Heitler-Matthews model is further modified fo the 
extended Heitler-Matthews model presented in Ref. [Q]. The extended Heitler-Matthews model 
includes the better description of pion interaction lengths and branching multiplicities. However, 
all of the above mentioned analytical models do not properly describe the shape of neither the 
electromagnetic nor the hadronic component. The shape of the electromagnetic component can 
be well reproduced by the intermediate shower model that modifies the original Heitler model, 
as explained in Ref. |^. In this study, we perform a similar alteration for the extended Heitler- 
Matthews model. Because of the complexity of the hadronic component we use a semi-analytical 
approach based on the branching probabilities and a simple numerical integration. 

2. Method 

For the description of the hadronic component we adopt the extended Heitler-Matthews model, 
see Ref. 0. A basic idea of the extended Heitler-Matthews model is visualized in Fig.|^. A pri¬ 
mary particle interacts in the atmosphere and produces secondary charged and neutral pions. The 
hadronic component within the model consists only of pions. Neutral pions immediately split into 
two photons and produce the electromagnetic showers that we do not examine. Charged pions fur¬ 
ther interact in the atmosphere and create other pions. A multiplicity Mji of produced pions in the 
pion initiated interaction in the shower is energy dependent [Q] 

( 2 . 1 ) 

E is the energy of the initial pion. We assume the ratio between charged and neutral pions to be 
2 : 1 in all interactions. Because of that the multiplicity of charged pions is related to the total pion 
multiplicity by 

( 2 . 2 ) 
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Figure 1: The hadronic cascade for energy dependent interaction lengths in the extended Heitler-Matthews 
model. Figure is taken from Ref. [Q]. 


The pion multiplicity of the first interaction depends on the mass of the primary particle. 
It differs for primary protons and iron nuclei. The pion multiplicity of the proton interaction is 
assumed to be equal to the pion multiplicity of the pion initiated interaction in the shower 


Mp«M: 

, = 0.15-£o°‘^. 

(2.3) 

For an iron primary it is considered to be 



Mpe - 

s0.45-£'^*^ 

(2.4) 

where Eq is the energy of the primary cosmic ray particle. 


The most important simplification of the model is the presumption that the energy is equally 

divided among all produced particles in the interaction. Therefore the pion energy after n interac- 

tions (energy branching) is given by 



F — 

^0 

(2.5) 


Mp/Pe-Mri’ 

Other necessary parameters of the model are 
nuclei we adopt [Q] 

the interaction lengths. 

For pions, protons and iron 

A;r-air[gCm^^]. ! 

=5^ 200-3.31og£'[eV], 

(2.6) 

Ap_air[gCm^2] : 

145-2.31og£'o[eV], 

(2.7) 

Ape-airfecm^^] : 

12. 

(2.8) 


In the extended Heitler-Matthews model 0, pions in a shower produce a cascade until they 
reach a critical energy and then decay. We do not use this assumption. We presume that interaction 
and decay processes of pions compete. To this end, we assign probabilities for the pion decay, for 
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their interaction with the air nuclei or their survival. When they interact they produce Mji pions 
(Eq.^jj) with energies of (see Eq.^3|). These assumptions allow us to calculate the number of 
pion decays in different atmosphere depths. Muons are produced at the same atmospheric depths 
where pions decay. 

In calculations we use an approximate relation between the depth in the atmosphere X and 
height h 


X(/i)[gcm 


1030•exp 


/r[km] \ 


(2.9) 


We model the atmosphere as N slant depths separated by a slant depth AX. We start from N = Q 
where the primary particle interacts and calculate the number of secondary pions and their ener¬ 
gies. Then we estimate the probabilities that pions decay, interact or survive until they reach the 
N = \ atmosphere depth. We assume that corresponding fractions of pions undergo the processes. 
Repeating all the calculations for the next N, and so on until the ground, gives us the energy dis¬ 
tribution of pions in the atmosphere. Only several energy bins are present because the energy is 
equally divided in interactions. Since the interaction and decay probabilities depend on the depth 
in the atmosphere^ and on the pion energy, we calculate only several probability values for each 
atmospheric depth. 

To reproduce the muon production depths (MPD) we have to incorporate the muon decay. In 
our model, as it is in relevant experiments [P, these production depths are determined for muons 
that reach the ground. Thus, for each muon, we calculate the probability that it arrives at the 
ground and assume this probability gives us the fraction of muons that contributes to the MPD. 
Energy losses of muons in the air are neglected. 


3. Results 

The comparison between the depths of muon production in the atmosphere and production 
depths of muons that reach the ground is depicted in Eig.p As shown, not only the total number 
of muons but also the shape of the profile and its maximum are affected by the muon decay in the 
atmosphere. We have to include the decay process to get reliable results. 

The distribution of energy branchings at which the pion decay contributes to the MPD is shown 
in Eig.p It is clear that for the example given in Eig.^ the majority of muons is generated at the 4* 
branching but the and 5* can not be neglected, otherwise the shape of the MPD profile will be 
affected. 

The most important results are MPD profiles calculated for different primaries. As an example, 
in Pig.@ MPD profiles for proton and iron primaries with energy of 10 EeV are presented. The 
maximum of the profile at the X^ax for the proton induced shower is larger and occurs deeper in the 
atmosphere than the maximum of the profile for the primary iron. Profiles in general correspond to 
those presented in Ref. [|]]. 

The energy evolution of the X^ax is shown in Pig.§. The difference between proton and iron 
showers is clearly visible and the general trend of the increasing X^ax as a function of energy 
is present. As a result of muon decays in the atmosphere, we can see unexpected features in 

*The zenith angle of the incident particle has to be taken into account. 
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Muon production profile, proton shower, E = 10 EeV, 0 = 60° 


xio^ 



Figure 2: Distribution of the production depth of muons. Green and blue curves represent respectively all 
produced muons and those that reach the ground. Results are for a primary proton with energy of 10 EeV 
and incident zenith angle of 60°. 


Contributions to MPD, proton shower, E = 10 EeV, 0 = 60° 



Figure 3: Contributions of pion decays at different energy branchings to the MPD for the example given in 
Fig.2. 
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Figure 4: MPD profiles for primary proton and iron particles with incident zenith angle of 60° at energy of 
10 EeV. 


this dependence. It is because muons are considered to have discrete energies at the production. 
Different energies contribute differently to the MPD as is visualized in Fig]^ The muon decay 
results in a suppression of lower energy branchings which leads to the decreasing of the X^ax when 
the higher energy branching starts to dominate the MPD. 


4. Conclusions 

We utilized the idea of the extended Heitler-Matthews model ii to calculate the muon pro¬ 
duction in the extended air showers. Despite the simplicity of the model we were able to reproduce 
MPD profiles for proton and iron primaries. Values of X^ax increase as a function of energy as 
expected but provide a systematic underestimation. The flattenings in the X^ax energy dependence 
show us the limits of the description and suggest that a more realistic energy division in the inter¬ 
actions is necessary for precise predictions of the X^ax- Details of the MPD calculations in real 
experiments are complicated and there is no straightforward way to directly compare them with 
our results. 
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Xmax energy dependence 



Figure 5: Energy dependence of the V^ax for protons and iron nuclei. Incident angles of the primaries are 
60°. 
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